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JE 要 本 文选 择 参数 化 的 电离 层 热 层 理论 模型 TIEGCM 作为 背景 模型 ， 基 于 
COSMIC 掩 星 观测 的 电子 密 ; 度 廓 线 数据 ， 应 用 集合 卡尔 曼 滤 波 方法 建立 了 全 球 电 
离 层 电子 密度 同化 模型 ， 实 现 了 全 球 电离 层 的 电子 密度 同化 。 同 化 结果 表明 ， 该 
E cc i.c 
与 背景 模式 相 比 ， 同 化 得 到 的 电子 密度 相对 于 观测 值 的 偏差 显著 下 降 。 此 外 ， 
组 同化 与 同时 同化 的 结果 对 比 显 示 , 平均 偏差 改善 基本 Paca ied 
偏差 在 峰值 高 度 以 上 略 有 减 小 。 

关键 词 电离 层 ， 数 据 同 化 ， 集 合 卡 尔 曼 滤波 ， 同 时 同化 
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Abstract Parameterized ionosphere model TIEGCM is used as background model. 
Basing on the COSMIC observations, global ionospheric electron density assimilation 
model is established using ensemble Kalman filter. Result shows, this model can 
effectively assimilate observations into background model and acquire three 
dimensional ionospheric electron density. Compared to background, the error between 
analysis and observations decreases significantly. The root mean square 
errorrRMSE) of NmF2 decreases about 60% for observations assimilated, and 20% 
for observations not assimilated. The RMSE of hmF2 doesn't get improvement except 
for mean error. 


The results of simultaneous assimilation (SA) and batches assimilation (BA) are 
compared for this case. The time that the two methods spend in assimilation is about 6 
to 7 minutes, which does not differ very much. SA needs nearly 8 GB storage while 
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BA less than 2GB. The statistic of electron density error shows that they nearly 
acquire the same mean error, but the SA gets relative better improvement in RMSE 


above 250km. 


Key words ionosphere, data assimilation, ensemble Kalman filter, simultaneous 


assimilation 


1 引言 


电离 层 的 电 


内 部 和 外 部 驱动 


化 而 变化 目 。 由 于 观 涡 


子 密度 分 布 可 以 影响 导航 、 短 波 通讯 系统 等 ， 
因素 影响 ， 随 经 度 、 纬 度 、 高 度 、 地 方 时 、 季 节 、 太 阳 活 动 及 地 人 厂 活 动 的 变 
上 缺乏 足够 的 时 空 分 辨 率 ， 仅 通过 观 涡 


电离 层 模式 
经 验 模式 基于 大 


1 难以 监测 电离 层 状态 。 


以 致 影响 人 类 活动 。 电 离 层 受 


直 以 来 是 进行 电离 层 研究 和 预报 的 基本 工具 ， 包 括 经 验 模式 和 理论 模式 。 


准确 的 短期 预报 
不 可 能 全 面 准确 


部 分 经 验 模式 和 
估计 ， 不 能 准确 


。 理 论 模式 基于 物理 方程 进行 建 模 ， 可 以 进行 短期 预报 ， 但 
的 考虑 电离 层 里 面 的 物理 过 程 及 扰动 情况 , 再 加 上 边界 条 件 和 初始 
确定 性 ， 及 数值 方法 带 来 的 误差 ， 理 论 模 式 进 行 的 预 


量 的 观测 数据 进行 统计 分 析 得 来 , 可 以 较 准 确 的 进行 中 长 期 预报 , 很 难 进行 


是 由 于 到 


理论 模式 可 以 模拟 电离 层 气候 特征 ,但 由 于 其 缺乏 对 电离 层 驱 动因 


的 进行 电离 层 现 报 、 预 报 品 。 


为 了 满足 对 


同化 技术 ,被 人 们 引 月 
多 种 电离 层 观测 ， 


x RO EO, 


地 基 GNSS-TEC 观测 和 天 其 掩 星 观 测 的 空 


漠 、 极 区 等 其 他 


观测 较 少 的 区 域 提 供 三 维 电子 密度 分 布 。 


国际 上 电离 


实验 室 一 南 加 州 大 学 全 球 同化 电离 层 模式 (JPL/USC GAIM ) 和 美国 
观测 同化 模式 (USU GAIMD P^. usu GAIM 的 同化 模式 又 包括 Gauss-Markov Kalman 滤波 同 
态 Kalman 滤波 同化 模式 站。 近年 来 国外 一 些 学 者 在 电离 层 同化 、 电 离 层 驱 


化 模式 和 简化 状 


层 数据 同化 研究 从 上 个 世纪 就 已 经 开始 , 其 中 比较 著名 的 , 如 美 
犹他 州立 大 学 全 球 电离 层 


B, dT AUS BU 


E 论 模式 
条 件 的 不 
常常 与 实际 情况 有 很 大 的 差别 中。 大 
素 的 可 靠 


电离 层 进行 监测 和 预报 的 需求 , 在 气象 上 已 经 运用 成 熟 的 技术 手段 一 一 数据 
昌 到 电离 层 研 究 和 预报 领域 。 在 过 去 的 十 年 上 
许多 电离 层 数据 同化 模型 被 开发 出 来 。 由 于 这 些 模型 主要 同化 有 限 空 间 履 
因此 并 不 能 精确 地 表征 全 球 范围 内 的 电子 密度 分 布 。 但 随 着 GNSS 的 发 展 ， 
间 履 盖 越 来 越 广 。 尤 其 是 掩 星 观 测 ， 可 为 海洋 、 沙 


LE 论 模式 和 


国 喷气 动力 


ES 


动 参数 以 及 电离 


层 热 层 的 相互 作用 对 电离 层 同 化 的 影 


了 不 少 成 果 。 国 内 这 方面 的 研究 还 不 是 很 多 , 主要 是 乐 新 安博 


拟 和 数据 同化 


究 , 其 在 电离 层 电 子 密度 空间 相关 性 方面 的 工作 在 电离 
大 的 应 用 价值 中。 
本 文 将 采用 集合 卡尔 曼 滤 波 同 化 方法 ， 把 COSMIC 电子 密度 廊 线 同化 到 NCAR 开发 的 


电离 层 
Hg 


和 等 方面 也 开展 了 一 些 研究 & D, pug 
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TIEGCM 模型 中 ， 给 出 初步 结果 ， 并 作 进 了 一 步 讨 论 ， 同 时 对 两 种 同化 方式 进行 了 比较 。 


2 同化 模型 构建 方法 


2. 1 集合 卡尔 曼 滤 波 方法 


Evensen 等 提出 基于 蒙特 卡 洛 随机 抽样 的 卡尔 曼 滤 波 方法 ， 也 就 是 集合 卡尔 曼 滤 波 


(Ensemble Kalman Filter, EnKF) 。 相 对 于 变 分 方法 和 


报 的 有 限 样本 来 


模 


非 线 性 


EAKS& JEJE, EnKF XJ 


EA 


A 


N O 


Ti 


估算 误差 协 方差 , 它 可 以 应 用 在 高 度 非 线性 的 动力 模型 中 ,， 既 不 需要 切线 性 


算 子 或 伴随 方程 ， 也 不 需要 随时 间 向 后 积分 。 其 基本 滤波 公式 如 下 : 


"i 
p ora (1) 

N—1 
K = LH'(HP,H! + R) 1 (2) 
X? = X? + K(Y — HX”) (3) 


式 中 X' 为 背景 扰动 场 ; P. EROR 22S Re 是 观测 误差 协 方差 矩阵 。K 为 增益 矩阵 ， 
用 来 修正 观测 点 附近 的 背景 场 。X&2 为 背景 场 ; 了 为 观测 向 量 ; 是 观测 算 子 ; Xe 为 同化 结果 ; 
EnKF 使 用 样本 统计 计算 背景 误差 协 方差 ， 容 易 产 生 虚 假 的 伪 相 关 ， 因 此 需要 用 相关 距离 来 
控制 伪 相 关 ， 即 协 方差 局 地 化 。 
2. 2 同化 变量 的 设置 

由 (3) 式 可 知 ， 同 化 所 需 变量 有 背景 场 X2、 观 测 值 Y、 由 背景 场 映 射 到 观测 值 的 插值 
算 子 或 观测 算 子 五 以 及 背景 场 、 观 测 值 的 误差 协 方差 矩阵 已 R。 下 面 介绍 这 几 个 量 的 获取 方 
法 。 


背景 场 由 TIEGCM 模式 预报 得 到 。NCAR 开发 的 TIEGCMI 是 一 种 三 维 、 时 变 、 地 球 热 层 
电离 层 耦 合 的 理论 模型 。 该 模型 使 用 有 限 差 分 技术 , 自治 求解 三 维 中 性 成 分 和 离子 的 动量 方 
程 、 能 量 方 程 和 连续 性 方程 。 其 默认 设置 的 经 度 、 纬 度 方 向 上 的 分 辨 率 分 别 是 5 度 和 5 度 ， 
高 度 方向 有 29 层 ， 间 隔 为 半 个 标高 ， 下 边界 约 为 97km， 上 边界 随 太阳 活动 不 同 在 500km 
到 700km 之 间 变 化 。 

TIEGCM 模式 是 一 种 电离 层 热 层 电动 力学 大 气 环 流 模型 ， 其 包含 很 多 大 气 参 量 的 预报 ， 
如 电子 密度 、 大 气温 度 等 。 这 里 选取 2006 年 第 80 天 的 数据 作为 初 值 ， 通 过 模式 预报 到 第 
159 天 0 时 0 分 。 然后 扰动 模式 输入 参数 的 太阳 活动 性 指数 F10.7, 构造 均值 为 150、 标 准 
偏差 为 50% 的 100 个 样本 集合 。 这 100 个 样本 作为 输入 参数 输入 到 模式 中 , 模式 预报 6 个 小 
时 ， 得 到 背景 场 X?。 

观测 值 为 cosMiC 掩 星 电 子 密度 廓 线 数据 。 这 里 所 用 的 COSMIC 掩 星 观 测 的 电子 密度 

亡 线 数据 是 UCAR 提供 的 二 级 电离 层 数据 产品 。 研究 表明 , 由 于 Abel 反 演 中 的 球 对 称 假设 ， 
F 区 高 度 以 下 的 电子 密度 的 反 演 结果 应 当 谨 慎 使 用 四。 所 以 本 文 只 对 150km 以 上 的 电子 密度 
进行 同化 ， 插 值 网 格 间距 为 10km. HX 2007 年 第 159 天 6 点 前 后 半 小 时 的 观测 数据 作为 观 
测 场 ， 对 观测 数据 的 质量 控制 ， 取 观测 值 大 于 零 的 观测 。 经 初步 质量 控制 后 约 有 3580 个 观 
测 点 ， 共 121 KAR. 

TIEGCM 输出 的 网 格 高 度 是 等 位 势 高 度 ， 为 了 便于 构造 观测 算 子 ， 这 里 将 位 势 高 度 换 
算 为 几何 高 度 , 并 插值 到 等 几何 高 度 面 。 观 测算 子 由 观测 点 所 在 网 格 立 方 体 空间 线性 插值 得 
到 。 

背景 场 的 误差 协 方差 由 集合 样本 统计 得 到 。 在 大 部 分 同化 系统 中 , 观测 误差 假设 是 相互 
独立 的 5 切 ， 本 文中 也 采取 同样 的 假设 ， 对 观测 值 不 考虑 协 方差 ， 假 设 观 测 偏 差 为 观测 值 
的 20%， 则 观测 误差 协 方差 如 下 所 示 : 


2 n de peg 
0, izj 


(4) 


2.3 局 地 化 

在 集合 卡尔 曼 滤波 方法 中 ,， 有 限 集合 数 造成 不 满 秩 和 估计 的 相关 虚假 噪音 。 为 了 消除 距 
观测 较 远 的 伪 相 关 ， 需 采取 协 方差 局 地 化 。 一 般 假 设 电离 层 空间 相关 性 满足 高 斯 分 布 ， 但 相 
关 距 离 目前 并 没有 很 好 的 结论 吗 。 电 离 层 的 相关 距离 可 以 看 作 是 各 向 异性 的 。 这 里 把 电离 
层 的 相关 性 简单 的 分 解 为 子午 方向 、 纬 线圈 方向 和 高 度 方向 的 乘积 中 。 电 离 层 的 相关 尺度 在 


m 


不 同 的 方向 有 显著 不 同 ， 且 有 明显 的 纬度 、 高 度 、 地 方 时 、 季 节 、 太 阳 活 动 变化 。 垂 直方 向 
的 相关 距离 随 高 度 升 高 而 增 大 ， 且 向 上 的 相关 距离 大 于 向 下 的 相关 距离 ， 总体 上 来 说 ， 纬 向 
相关 距离 比 子午 向 相关 距离 大 。 

把 集合 样本 估计 的 背景 误差 协 方差 和 局 地 相关 函数 相 乘 ， 得 到 新 的 协 方差 矩阵 。 


pij = e *XGi* /L*) x e ex (Gu /Ly’) (5) 


Pij 为 6j 两 点 的 相关 系数 ，qi、0ij 分 别 表示 两 个 点 在 经 度 和 纬度 方向 的 实际 距离 。 若 定义 两 

个 点 之 间 的 相关 系数 为 0.75 时 的 距离 为 相关 距离 , 根据 乐 新 安 等 外 人 的 研究 , 纬度 方向 的 相 

关 距 离 Lx 取 10”， 经 度 方向 的 相关 距离 Ly 从 中 纬 ~40” 变 化 到 赤道 地 区 ~20”， 高 度 方向 设 

置 不 相关 。a 为 相关 距离 系数 ，0.75 = e-%,a 为 0.2877. 

2.4 同化 方式 

顺序 同化 是 对 观测 数据 逐个 进行 同化 , 同时 同化 是 同时 对 所 有 观测 数据 进行 同化 。 顺序 

l 同化 是 针对 观测 数据 的 局 部 优化 , 同时 同化 是 针对 所 有 观测 数据 的 全 局 优化 。 ETE RU, 
> 当 观 测 误差 不 相关 时 , 在 集合 卡尔 曼 滤波 中 顺序 同化 和 同时 同化 差异 不 大 。 顺序 同化 每 次 只 
= 同化 一 个 观测 数据 , 对 计算 机 的 内 存 大 小 要 求 低 , 但 将 所 有 观测 数据 都 同化 完 , 所 需 时 间 长 。 
同时 同化 一 次 同化 所 有 观测 数据 ， 计 算 效 率 高 , 但 要 求 计算 机 的 内 存 大 , 被 同化 的 数据 量 


时 ， 内 存 大 小 会 无 法 满足 计算 需求 。 将 顺序 同化 和 同时 则 化 的 优点 相 结合 ， 对 观测 数据 进行 
分 批 同化 ， 即 为 分 组 同化 。 
co 为 比较 分 组 同化 与 同时 同化 的 差异 ,本 文 分 别 用 这 两 种 方法 进行 同化 , 并 对 同化 结果 进 
> 行 了 对 比分 析 。 本 文 进行 分 组 同化 时 ， 采 取 逐 个 观测 廊 线 成 批 同化 的 方法 。 


` 3.1 同化 结果 
在 121 条 观测 廓 线 中 随机 取出 50% 进 行 同化 , 同化 方式 为 同时 同化 , 利用 剩余 50% 的 未 
参与 同化 的 廓 线 来 验证 同化 效果 ， 下 面 给 出 了 一 些 结果 。 
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图 1 同化 前 后 背景 
实 线 为 同化 结果 。 


场 、 分 析 场 与 参与 同化 的 电子 密度 廓 线 的 比较 。 点 画 


线 为 观测 场 、 虚 线 为 背景 场 、 


Fig.1Comparison of electron density profile between background (dash lines) and EnKF assimilated (solid 


lines) for assimilated obs 


ervations (dash-dotted lines). 
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2 同化 前 后 背景 场 、 


分 析 场 与 未 参与 同化 的 电子 密 / 


Fig.2Comparison of electron density profile between background (dash lines) and EnKF assimilated (solid 


度 廊 线 的 比较 。 


lines) for unassimilated observations (dash-dotted lines). 


这 里 给 


浅 灰 色 虚 线 为 背景 场 、 黑 色 实 线 为 分 析 场 。 
果 。 从 同化 结果 可 以 看 出 ， 
未 参与 同化 , 但 分 析 场 仍 与 观测 场 符合 得 更 好 , 这 
的 同化 效果 , 说 明 该 同化 模型 能 够 将 观测 资料 有 效 地 同化 到 背 
密度 分 析 场 ， 对 没有 观测 资料 的 地 方 也 能 


测 点 处 的 结 
DUIS 


3. 2 相对 偏差 


出 了 同化 结果 与 观测 场 及 同化 之 前 的 背景 场 的 比较 个 例 , 深 灰 色 点 画 线 为 观测 场 、 
参与 比较 的 背景 场 和 分 析 场 廓 线 为 插值 到 相应 观 


有 很 好 的 改善 作用 。 


图 1 中 分 析 场 与 观测 值 更 为 接近 ,图 2 
文 是 由 背景 场 提供 的 相关 性 所 带 来 
9 景 模式 中 ,获得 全 球 三 维 电子 


中 虽然 这 些 


为 进一步 验证 同化 效果 , 统计 了 同化 前 后 在 不 同 高 度 上 背景 场 和 
差 。 相 对 偏差 、 平 均 偏 差 和 和 标准 偏差 公式 分 别 为 : 
e; = 5 x 100% 
Yi 
e 三 站 ei 


Xi- (ei-8)? 


O. = 
xX n-1 


同化 结果 相对 观测 的 偏 
(7) 


(8) 


(9) 


中 


差 ， 


ei 为 电子 密度 的 相对 偏 


e 和 ax 分 别 为 对 应 高 度 上 电子 密度 相对 1 


中 ，X%i 为 某 位 置 上 的 背景 场 或 同化 结果 ， 区 为 对 应 位 置 的 观测 值 ， 


局 差 的 平均 偏差 和 标准 偏差 。 
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Fig.3.The statistic of electron density error at observation profiles assimilated for simultaneous assimilation. 


(Left: before assimilation, right: after assimilation). 
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为 同 


兰 统 计 《〈 左 图 为 同化 前 ， 右 


ET ER 


4 同时 同化 前 后 未 参与 同化 观测 廓 线 处 的 各 个 高 度 上 
化 后 ) 


Fig. 4.The statistic of electron density error at observation profiles unassimilated for simultaneous 
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assimilation. (Left: before assimilation, right: after assimilation). 
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场 样本 的 标准 偏差 相对 于 样本 均值 的 百分比 


5 东经 115 度 、 北 纬 40 度 处 背景 
Fig. 5 The profile of the relative standard deviation of background on 115°E 40°N. 


局 差 ， 两 条 虚线 分 别 为 平均 侦 差 加 减 一 


D 


图 3、 图 4 左 图 中 的 实 线 为 同化 前 背景 场 的 平均 1 
倍 的 标准 偏差 。 右 图 为 同化 后 分 析 场 的 偏差 统计 。 可 以 看 出 无 论 是 参与 同化 还 是 未 参与 同化 
前 背景 场 的 平均 偏差 和 标准 偏差 均 有 所 减 小 ， 


的 观测 廓 线 处 , 总 体 上 同化 后 分 析 场 相 比 同化 之 
在 300km 以 上 平均 偏差 由 同化 前 的 50% 降 低 到 同化 后 的 15% 以 内 ,标准 偏差 由 同化 前 的 50% 


图 中 也 发 现 ， 在 较 低 高 度 同化 效果 并 不 理想 ， 


一 是 由 于 


以 上 降低 到 同化 后 的 40% 以 内 。 从 两 
此 高 度 范围 内 的 背景 值 偏 小 , 男 一 个 原因 是 背景 场 的 样本 集合 在 此 范围 内 样本 偏差 较 小 造成 
的 。 图 5 给 出 的 是 东经 115 度 、 北 纬 40 度 处 背景 误差 协 方差 对 角 线 值 的 平方 根 相 对 于 集合 
样本 平均 值 的 百分比 ， 从 图 中 可 看 出 ，250km 以 下 由 集合 样本 统计 获得 的 背景 误差 在 2096 
以 下 ， 远 小 于 350km 以 上 的 背景 误差 。 背 景 值 与 相对 偏差 偏 小 导致 背景 误差 偏 小 ， 同 化 结 
果 就 会 更 接近 背景 场 ， 与 观测 数据 之 间 的 偏差 改善 就 会 大 大 减 小 。 
3. 3 相关 性 比较 
5 NmF2 5 NmF2 
15 x10 15 x10 
corrcoef 0.88 corrcoef 0.98 
mean error 37.496 mean error -0.1% 
o std error 87.8% © (7 std error 26.8% © 
E 10 © 5 [e 
o o 
z 三 
0 Qo 
[9] [9] 
Ó Ó 
10 15 10 15 
backgroud (cm?) x10° analysis (cm ?) x10 
hmF2 hmF2 
400r corrcoef 0.79 NH > 400r corrcoef 0.69 P 
mean error 44.1km + mean error 24.5km kA 
EM + E -Z 
€ 300} std error 28.2km E L} £30 std error 37km F gi 
Pa imd Ar T 
2 KG 9 + 性 T g 
& 200: ^ = 200; pi 
O 〇 OO 
OO Ó 
100 100 
0 "d " A É é 0 一 r r r r 
0 100 200 300 400 0 100 200 300 400 
background (km) analysis (km) 


D 


6 同化 前 后 参与 同化 观测 廓 线 处 NmF2 和 hmF2 的 相关 性 统计 


Fig.6 Correlation statistic in NmF2 and hmF2 (with assimilation). 
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图 7 同化 前 后 未 参与 同化 观测 廊 线 处 NmF2 和 hmF2 的 相关 性 统计 
Fig.7Correlation statistic in NmF2 and hmF2 (with no assimilation). 
图 6、 图 7 中 分 别 给 出 了 电离 层 峰 值 电子 密度 NmF2 同化 前 的 背景 场 与 观测 值 、 同 化 后 
分 析 场 与 观测 值 的 散 点 图 。 左 图 横 坐 标 为 背景 场 峰值 电子 密度 ， 纵 坐标 为 观测 值 ; 右 图 横 坐 
图 中 可 以 看 出 对 于 参与 同化 的 观测 ， 同 化 前 的 背景 场 与 


标 为 分 析 场 ， 纵 坐标 为 观测 数据 。 从 
观测 值 的 相关 系数 为 0.88， 同 化 后 为 0.98， 同 


的 37.4% 下 


偏差 并 未 改善 。 
3.4 同时 同化 与 分 组 同化 


攻 后 相关 系数 有 


>> 


E 想 ， 平 均 偏差 有 所 下 降 , 但 相关 系数 和 和 


民 大 提高 , 平均 偏差 从 同化 前 
降 到 -0.1%， 标 准 偏差 降低 约 60%。 对 于 未 参与 同化 的 观测 ， 相 关系 数 从 同化 前 的 
0.81 提高 到 0.94, 平均 偏差 从 同化 前 的 17.4% 减 小 到 3.6%， 标准 偏差 减 小 约 20%。 从 图 


值 高 度 比 较 结果 来 看 ， 峰 值 高 度 同化 结果 并 不 型 


中 峰 
标准 


= Mean error 
=== Mean error + Standard error 
=== Mean error - Standard error 


= Mean error 
==- Mean error + Standard error 
==- Mean error - Standard error 


450} iy? 450} Fi 
1 1 Y 
1 1 1 
/ Y M 
400 ' i 400 ( 
I 1 EY 
[ H 1 
l A m 
— 350} 全 本 一 350| / 
E I 1 E / 
r3 M E z i 
+= 1 += 
E 300 Iji! E 300 ( 
D 1 1 o V 
I / li I 1 
IN 1 y I 
250} /! (d. 250} t 
1 lil £ 
Á 1 1 t 
2 1 1 I 
200 1 1 1 200 t 
) 
" 17 
} Ii | i 
-150 -100 - 50 100 150 -150 -100 - 50 100 150 
Relative error of electron density (%) Relative error of electron density (%) 


pm 图 8 同时 间 化 后 各 个 高 度 上 电子 密度 误差 统计 与 分 组 同化 对 比 。( 左 图 为 同时 同化 , 右 图 为 分 组 同化 ) 


Fig.8 Comparison of the statistic of electron density error between different assimilation methods; the left is 


simultaneous assimilation and right batches assimilation. 


将 所 有 观测 数据 同时 同化 到 背景 场 中 需要 很 大 的 内 存 存储 , 该 试验 是 在 多 核 计算 机 上 进 
Co 行 的 ，CPU 主 频 为 1.87GHz， 同 化 的 观测 点 数 为 3579。 同 时 同化 最 大 占用 内 存 约 8GB， 同 化 
N 耗 时 6 分 36 秒 。 为 比较 同时 同化 与 分 组 同化 的 差异 ， 对 同样 的 观测 数据 采用 分 组 同化 进行 
m 了 同化 实验 ， 而 分 组 同化 占用 内 存 268 以 内 ， 耗 时 7 分 21 秒 。 如 图 8 所 示 ， 同 化 结果 对 比 
= 发 现 ， 同 时 同化 和 分 组 同化 相对 观测 值 的 平均 偏差 基本 一 致 ， 但 同时 同化 后 的 标准 偏差 在 
ku 300km 高 度 以 上 相对 更 小 。 由 此 可 见 ， 同 时 同化 具有 全 局 优化 的 优点 , 但 占用 内 存 大 ; 分 组 
一 同化 能 够 取得 不 错 的 同化 效果 ， 且 占用 内 存 小 。 


4 结论 与 展望 

= 本 文 给 出 了 一 个 基于 集合 卡尔 曼 滤 波 技 术 的 电离 层 数据 同化 模型 的 初步 构建 和 试验 结 
€ 果 。 我 们 利用 集合 卡尔 曼 滤 波 同化 方法 , 将 COSMIC Ti ELI BU Ha T A E E 2x [8] £558] TIEGCM 
背景 模型 中 。 同 化 结果 表明 ， 采 用 集合 卡尔 曼 滤波 算法 ， 把 电子 密度 同化 到 背景 模式 中 能 够 
获得 较 好 的 同化 结果 ， 对 于 参与 和 未 参与 同化 的 观测 ，NmF2 的 标准 偏差 分 别 降 低 约 60961 
20%， 说 明 我 们 设计 的 算法 可 行 、 所 选择 的 各 种 参数 比较 合理 。 分 组 同化 与 同时 同化 的 对 比 
结果 表明 ， 分 组 同化 对 内 存 需 求 小 ,但 同时 同化 结果 更 好 ， 且 耗 时 更 短 ， 条 件 允 许 的 情况 下 
宜 采 用 同时 同化 的 方法 。 
如 前 所 述 ， 电 离 层 驱动 参数 对 电离 层 预 报 是 至 关 重 要 的 。 仅 同化 获取 电子 密度 ， 没 有 优 
化 的 电离 层 驱 动 参数 ， 难 以 对 电离 层 进行 可 靠 的 预报 。 本 文 构建 的 电离 层 同 化 模型 ， 只 可 用 
于 电离 层 现 报 , 接 下 来 我 们 将 进行 电离 层 驱 动 参数 的 同化 预报 研究 。 此 外 ， 电 离 层 同化 预报 
还 涉及 到 观测 数据 的 质量 控制 、 超 大 矩阵 的 存储 与 计算 、 误 差 协 方差 的 选择 等 问题 ， 还 需 对 
算法 进行 进一步 优化 ， 以 改善 同化 结果 的 精度 、 提 升 同 化 速度 。 
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